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We devise a phase-coherent three-pulse protocol to probe the polariton dynamics in a 
trapped-ion quantum simulation. In contrast to conventional nonlinear signals, the pre¬ 
sented scheme does not change the number of excitations in the system, allowing for the 
investigation of the dynamics within an A^-excitation manifold. In the particular case of 
a filling factor one (N excitations in an A^-ion chain), the proposed interaction induces 
coherent transitions between a delocalized phonon superfluid and a localized atomic insu¬ 
lator phase. Numerical simulations of a two-ion chain demonstrate that the resulting two- 
dimensional spectra allow for the unambiguous identiflcation of the distinct phases, and the 
two-dimensional lineshapes efficiently characterize the relevant decoherence mechanism. 
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I. INTRODUCTION 


Polaritons - hybrid quasiparticles originating from the strong coupling between light modes 
and matter excitations - represent an important research area, both for fundamental and practical 
reasons: Their bosonic character may lead to Bose-Einstein condensatiorP^l, and the larger co¬ 
herence length of the light modes may enhance the carrier mobility in organic semiconductorP. 
While originating primarily in quantum-optical settings such as Josephson junctions and arrays 
of coupled cavitieP^, polaritonic excitations are now also studied in molecular systemP. Alter¬ 
natively, trapped ions can be used to study the properties of strongly correlated systems under 
well-controlled conditions and with manageable, slower time scales than solid-state systemP^. 

Polaritonic systems can be modeled in ion trap experiments by coupling the electronic and vi¬ 
brational degrees of freedom of the ions confined by a harmonic potential. The atomic excitations 
are encoded into the ions’ electronic states, and the local vibrational modes of the ions, described 
by quantum mechanical harmonic oscillators, take on the role of the light modes. The Coulomb 
repulsion between the ions leads to couplings of the local phonons. Finally, a suitably chosen 
external laser field realizes Jaynes-Cummings type interaction^^ between each ion’s electronic 
state and its local phonon modeP^i^. The full system is then described by a Jaynes-Cummings- 
Hubbard modepl, in which the total number of electronic and vibrational excitations is conserved. 

For filling factor one (N excitations on N sites), the ground state varies, with the detuning 
of the external laser field, from a phonon superfluid phase, in which vibrational excitations are 
delocalized over the entire chain, to an atomic Mott insulator phase, in which electronic excitations 
reside on individual ions. While all excitation energy is stored in the vibrational degree of freedom, 
in the former case, it is fully absorbed by the electronic degree of freedom, in the latter. In the 
intermediate regime, the system exhibits a polaritonic phase of coupled atom-phonon excitations. 
When the filling factor is varied, the model allows for further phases, such as a polariton glasP. 

In this work, we devise an interaction scheme which induces coherent transitions between these 
two phases. This interaction will allow us to design pulse sequence^^^^^^ for multidimensional 
nonlinear spectroscopy of the polariton dynamics. In contrast to conventional spectroscopic tech¬ 
niques, our interaction does not change the number of excitations in the system, but rather switches 
between electronic or phononic character of the excitations. In a phase-coherent three-pulse se¬ 
quence, the system evolves through both the superfluid and the insulator quantum phases, as well 
as a coherent superposition of both. We will see that nonlinear spectroscopy is particularly well 
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suited to identify the relevant excitation processes, as well as the dominant decoherence mecha¬ 
nisms via the observed two-dimensional lineshapes. To ease scalability towards large system sizes, 
in which quantum phase transitions become most pronounced, we renounce of single-ion address¬ 
ability in our spectroscopic protocol, which, however, could be exploited to add spatial resolution 
to the obtained nonlinear spectraP^. 

II. BACKGROUND 

A. Dynamics and laser control of trapped ions 

The dynamics of polaritonic systems can be simulated by subjecting a collection of trapped ions 
to appropriate laser field d^^l^^ i taking into account both the electronic and motional states of all the 
ions. As mentioned before, the relevant electronic states of the trapped ions are modeled as non¬ 
interacting spin-1/2 system^, while the motion is determined by their common trap potential, 
which, to a good approximation, is described by a three-dimensional harmonic oscillatoi^. The 
analogy to polaritonic systems is achieved by identifying the phonons in the ion trap with the 
photons of coupled cavities. In the remainder of this article we refer to the actual physical degree 
of freedom, i.e., the phonons, rather than the mathematically equivalent photons in the case of an 
optical polaritonic system. 

Near their motional ground state, the interplay between the potential and the strong Coulomb 
repulsion of the ions determines the chain geometry. The equilibrium positions of the ions are 
determined by the ratio of the trap frequencies fy^Pz) along the different spatial directions. 
A linear chain arrangement is achieved when, e.g. ^ Uy ^ Uz, and allows for convenient 
laser addressing of individual ions, with typical inter-ion spacings of 5 — 15 /im and typical trap 
frequencies on the order of Px^i^y ~ 1 — 10MHz and i/z ^ 0-f — IMHz. In the following, we 
will consider such a linear chain of laser-cooled ions, which allows us to assume that the ions’ 
displacement from their equilibrium positions is much smaller than their separation. Employing a 
second-order expansiorP^in terms of the displacement along the a:-direction leads to the Hamilto¬ 
nian (we set ^ = 1) 

N N N 

Ho = + X] {a\ai + ci\a^ , (1) 

k=l k=l k,l=l 

k<l 

where a\ creates a local phonon along the a:-direction at the ion k, and is the spin-1/2 lad- 
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der operator at site k. The electronic states of all the ions are described by the same transition 
frequency u. Due to large spatial separation between adjacent ions, direct transfer of electronic 
excitations between the ions may be neglected. However, the Coulomb repulsion of the ions affects 
the motional degrees of freedom in Eq. ([T]): It is responsible for both, modulations of the local os¬ 
cillator frequencies cok = — ^he phonon-phonon couplings tki = z/a,/5/(2|ufc — ui\^), 

where (3 = is the ratio between the confining trap frequencies in axial (z) and radial (a:) 

direction, and ui the ions’ equilibrium positions in units of Iq = (e^/ 

An essential tool for the control of the quantum state of trapped ions are coherent laser 
interaction^ ^ The interaction Hamiltonian, generated by a continuous-wave laser, may be 
written as 

^^int = ^{^t + ^k) , ( 2 ) 

and describes the interaction of the k-th ion (at the position fk) with a laser beam of frequency col, 
wavevector kL, and Rabi frequency Q. For laser-cooled ions of low kinetic energy, the trap poten¬ 
tial is well approximated by a static harmonic potentiaP^, and we may express the ion position in 
terms of the harmonic ladder operators. For example, when kL coincides with one of the motional 
axes, e.g., the ^-direction, in the case of a single ion, we obtain kiVk = {kil^k = v{<^k + a\), with 
the bosonic ladder operators as defined above, and the Lamb-Dicke parameter 77 = k/\/2mv^, 
which depends on the frequency of the harmonic motion. In general, for arbitrary numbers 
of ions, tuning the frequency lol allows to engineer effective couplings between the electronic 
and vibrational degrees of freedom, as well as between the electronic states of different ion^^^EH 
To see this, it is important to notice that the electronic spectrum of the ions is modulated by vi¬ 
brational sidebands [as can be seen in Eq. ([T])], which originate in the ions’ motion around their 
equilibrium positions. In analogy to Stokes and anti-Stokes lines in Raman spectroscopjP^, these 
resonances are found when the laser detuning coincides with multiples of the frequency of the 
harmonic motion. 

Using tunable narrow-band lasers, specific sidebands can be resonantly selected. This effec¬ 
tively creates coupled electronic and vibrational excitations: Driving the first red (blue) sideband 
leads to laser-induced transitions which create an electronic excitation while a phonon is removed 
(created). Hence, the dynamics of, e.g., the red sideband is described by a Jaynes-Cummings 
Hamiltonian which induces couplings of the type ~ cr^afc -|- aka\. A more detailed treatment that 
includes the contribution of various harmonic modes is described in appendix [A} 


4 







B. The polaritonic Hamiltonian 


In this paper, we consider the scenario considered in refM To couple motional and electronic 
degrees of freedom, and thereby create polaritonic excitations, one irradiates the entire ion chain 
with a travelling-wave laser, oriented along the radial x-direction, whose frequency is chosen as 
ujl = uj — Px + ^, i.e., the laser is detuned from the radial red sideband transition by Apl 

In the Lamb-Dicke limit (c.f. appendix |^, and after performing the rotating wave approxi¬ 
mation, the dynamics of N ions is effectively described in an appropriate rotating frame by the 
HamiltoniarfAl 

N N N N 

H = '^uJkaiak+Y^ tki (alai + ajak^ + A'^a+a^ + gY^ (a+ aj) , (3) 

k=l k,l=l k=l k=l 

k<l 

where, as before, a\ creates a local phonon along the x-direction at site k, and is the spin-1/2 
ladder operator of ion k. The coupling between spin and phonon degrees of freedom is denoted 
with an effective Rabi frequency g. The detuning A from the red sideband transition determines 
the effective local spin’s energy splitting. The validity of the above Hamiltonian is limited to the 
case where all parameters /?, A, and g are smalP^l, for the following reasons: The ratio of squared 
trap frequencies /? must be small to justify the quadratic expansion of the Coulomb interaction 
term which lead to the motional Hamiltonian ([T]). The detuning A must be significantly smaller 
than the trap frequency to ensure the validity of the rotating wave approximation which was 
employed to obtain Eq. 0. Finally, when the laser-induced atom-phonon coupling becomes too 
large, multiple phonons can be created or destroyed in one electronic transition, thus higher-order 
terms would have to be included in the interaction term in Eq. Q. 

We reiterate that the coupling between spins and phonons is of polaritonic (rather than pola- 
ronic) nature, since it leads to an exchange of atomic and phononic excitations: One easily verifies 
that the total number of excitations N = + cr/'cr/') is a constant of motion of the Hamil¬ 

tonian (|3]).^Due to the interaction between phononic and electronic degrees of freedom, an initial 
spin excitation localized on a specific ion will be transferred to the phononic subspace, spread due 
to the Coulomb interaction, and transfer back to the electronic subspace. Our present analysis 
aims to understand the dynamics in a given A-excitation manifold of the system. In particular, 
we restrict to the case where the number of excitations equals the number of ions, i.e., for a filling 
factor one. 

In this case, the ground state of the system’s A-excitation manifold for —A/g ^ 1 is given by 
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FIG. 1. (a) Schematic representation of the trapped-ion quantum simulation of polaritonic systems. Laser 
fields couple the electronic states of the ions to their motional degrees fo freedom, which in turn are coupled 
by the common trap potential, (b) Two-excitation manifold of the Hamiltonian Q for = 2 ions as a 
function of the ratio between detuning and interaction strength A/^. The eigenstates are indicated, their 
subscripts denote the number of spin excitations, and the superscripts label the individual states within a 
spin excitation submanifold. 

the uncorrelated atomic insulator (atl) state 

|atl) = I t)®^ ® |0)®^ (4) 

with localized spin excitations. Here, | t) denotes the excited electronic state, and |0) the phononic 
vacuum. 

Conversely, for A/^r ^ 1, the lowest-energy state is characterized by excitations of the mo¬ 
tional eigenmode of lowest energ)!^, created by the operator h\ = This state is denoted 

the phonon superfluid (phSF) state, 

IphSF) = I 4.)®" ® ^ (Ml)" 10)®", (5) 

since the excitations are spatially delocalized within the phonon degrees of freedom. In the inter¬ 
mediate regime. A/p 0, the system is described by a polaritonic superfluid phaseP^. This phase 
is characterized by a delocalized ground state of polaritonic character, i.e., containing both atomic 
and phononic excitation^. 

Let us now also briefly discuss the excited states for the simplest nontrivial example: the two- 
excitation subspace in a two-ion chain. The structure of the eight-dimensional spectral manifold 
is depicted in Fig.jT]?), as a function of the ratio A/p. 
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If A/^" ^ 1, the spin Hamiltonian A dominates over all other terms in the Hamil¬ 

tonian. In the case of the two-ion chain considered here, this Hamiltonian has three eigenvalues 
0, A, and 2 A (which are determined by the number of excited spins), and consequently we iden¬ 
tify three submanifolds of states, which differ in the number of spin excitations: The submanifold 
(l^^o )) 1 ^^)) 1 ^?)}’ with the ground state )’ encompasses all three states with close to zero 
spin excitations. The first excited-state submanifold contains four states |^f), 

with a single spin excitation. Finally, both spins are excited in the highest excited state |^^ 2 )- Since 
our analysis restricts to the two-excitation subspace, the respective spin excitations are comple¬ 
mented by phonon excitations among the two eigenmodes to reach a total of two excitations. 

In the other limit, when —A/g ^ 1, the same three submanifolds can again be identified, 
however, their energies are reversed: The state |\& 2 ) now constitutes the ground state manifold, 
and the zero-spin submanifold yields the highest-excited states. 

In between these two limits, when \A/g\ < 1, the eigenstates lose their clear character, and in 
the thermodynamic limit, one expects an abundance of avoided crossings which accompany the 
quantum phase transition and point to quantum chaotic behavioi^. 


C. Phase-coherent two-dimensional spectroscopy 

Nonlinear spectroscopy provides powerful tools to probe the dynamics and decay processes 
of complex interacting quantum system^. A general theoretical framework in the context of 
trapped-ion systems was recently developecP^, which allows for the identification of a variety 
of phenomena ranging from decoherence-free subspace^to phonon currents in a nonequilibrium 
steady- stateP^ and structural phase transitions in ion crystal^. The two essential ingredients are 
sequences of short, perturbative pulses separated by tunable time-delays, and a mechanism that 
allows to post-select the phase signature of particular excitation sequences. In large ensembles, 
the latter can be achieved via phase matching, which is not possible for trapped-ion systems, due 
to the limited amount of scattered photons per incoming probing puls^. A feasible alternative 
is presented by phase-cycling which was developed to probe nuclear magnetic resonance systems 
with nonlinear spectroscop}^^. 

In the following, we shall employ these ingredients to investigate the dynamics induced by 
Eq. 0. 
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III. RESULTS 


A. Probing the A^-excitation subspace with a phase-coherent three-pulse scheme 

For the present simulation study, we consider the interaction with a global laser field, which 
illuminates all the ions in the trap with equal intensity and is tuned resonantly to the motional 
sideband of the collective mode bi. This allows to realize an interaction of the forrrP^^ (see 
appendix [A|) 

V, = (6) 

= i(v,(4.)-Vl{4,)'j, (7) 

where J± = creates collective excitations among the spins, where the coefficients 

Cfci decompose the collective mode bi into local contributions as bi = Ckiak- It is possible to 
switch between the atomic insulator and the phonon superfluid states by inducing N interaction 
events of the above form, when bi describes the vibrational eigenmode of lowest energy. These 
interactions conserve the number of total excitations in the system: Each spin excitation is accom¬ 
panied by a phonon annihilation, and vice versa. Phase cycling with respect to the phase cj) can 
distinguish between the different contributions to the total signal. This technique further renders 
weak, perturbative pulses sufficient to obtain multidimensional spectra combined with a precise 
time resolution, and allows to post-select the evolution of specific coherences. 

Based on an extension of the single quantum coherence signaP^, we design the following exci¬ 
tation scheme (see Fig.|^: 

1. The system is initialized in the phonon superfluid state |phSF). 

2. During the first time interval, ti, the system evolves in the phonon superfluid state Q. 
This evolution is described by the Liouville space superoperator Q{t) = exp[£f], where the 
Liouvilian C contains the coherent evolution according to the Hamiltonian Q, as well as 
the influence of possible noise sources. Some of the dominant noise sources of trapped-ion 
experiments are mentioned in the course of our simulations in the next section. We obtain 
^(ti)|phSF)(phSF|. 

3. The first excitation pulse populates higher-excited states of the system. By phase cycling 

with respect to the phase signature — ((> 2 ), where (pi denotes the phase of the z-th 
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FIG. 2. Pulse sequence for the nonlinear measurement protocol: The system is prepared in the superfluid 
phase IphSF), and evolves freely during ti. Two excitation pulses then transfer the chain into the atomic 
insulator phase |atl). Post-selecting the phase N{(j )2 — 4>2) allows to separate the excitation process into 
two steps, and to study the evolution of the coherence |atI)(phSF| during t 2 . The readout signal is created 
at the end of the sequence by fluorescence detection at ion j. 

pulse [see Eq. (|^], it is possible to post-select those processes, in which the excitation pulse 
transfers the system into a superposition state of phonon superfluid and atomic insulator. 
Thus, the contributing terms are described by (vj^^)'^^(fi)|phSF)(phSF|, where = 
J+bip represents the Liouville space superoperators associated with an excitation acting 
only on the left (L) side of the density matrix. Similarly, we introduce Vj^^p — pJ-bl 
The decay of this superposition can be observed during the second time interval t 2 , and is 
described by g(fa) (vf (fi) |phSF) (phSF | 

4. The phase cycling procedure described above post-selects pathways, in which the second 
excitation transfers the system fully into a population of the atomic insulator state, whose 
dynamics is monitored during t^. We obtain 0(f3)(vj'^^)'^^(f2)(V}^^)^0(fi)|phSF)(phSF| 

5. Finally, the signal is obtained by measuring the spin population of one of the ions (charac¬ 
terized by the ion’s index j) via fluorescenceP^i^III] 

The measured signal is given by the nonlinear correlation functiorpl 

S{tl,t2,t3]j) = tr |cr+(Tj‘^(f3)(v}^^)^0(f2)(V}^^)^^(fi)po} , (8) 
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where po = |phSF)(phSF| denotes the initially prepared superfluid state and tr the trace. To 
efficiently analyze this multidimensional signal, we will Fourier transform two of the time delays, 


or 


/ dti 


= / dt2 




(9) 


( 10 ) 


Let us discuss the conditions and implications of each of the above steps in further detail. The ini¬ 
tialization of the phonon superfluid state can be achieved via a combination of ground-state laser 
cooling, optical pumping, and well-deflned pulses on blue sideband transition^^®21 jf system 
parameters are chosen such that the initial state |phSF) actually represents an eigenstate, there 
is no time evolution during ti. The presence of small contributions of other excited states will 
become visible in the resulting two-dimensional spectrum. Such contributions would also open up 
the possibility of removing spin excitations with the first pulse. For |A/p| ::§> 1, the state |phSF) 
is very close to an eigenstate of the system: When A > 0, it constitutes the ground state, when 
A < 0 the highest excited state, with small corrections due to the finite spin-phonon coupling g. In 
this case, we can neglect contributions from other diagrams, in which, for instance, an initial spin 
excitation creates additional pathways. These need to be taken into account only when \^/g\ < 1, 
which shall not be our concern in this paper. Finally, the measurement of spin populations is 
easily realized by fluorescence detectioiJAEMI] Since the spatial decomposition of the addressed 
eigenmode is contained in the spin excitation operator J, the signal can be dependent on j, i.e., 
we obtain spatial information about the relevant vibrational mode only with local readout, without 
resorting to local excitations. Single-site addressing is in principle possible and is able to enhance 
the information content compared to conventional nonlinear spectroscopic signal^^^^. Whereas 
the spatially resolved readout of the fluorescence light is routinely realized via high-resolution 
cameras, even for large N, the addressing of individual ions with focussed lasers becomes techni¬ 
cally demanding for increasing system size^. Since the macroscopic properties of the quantum 
phases and the corresponding phase transition is most interesting in this large-A limit, we avoid 
the usage of focussed excitations to ease scalability. 

In the present analysis, we have considered the impulsive limit, which implies that the length 
of the applied pulses is much shorter than the characteristic time scale of the relevant system 
dynamics. The dynamics of the effective polaritonic Hamiltonian Q is determined by the motional 
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FIG. 3. Center: Absolute value of the two-dimensional spectrum = 0, f22)f^3)| (where we omit 

the site index j in ( [T0| ), since the fluorescence signal is independent of j for = 2) of the pulse sequence ([^ 
in a two-ion chain implementing the Hamiltonian ([^, with Ux — ^ MHz, (3 — 5 kHz, A = 50 kHz, g — 
5 kHz. The panels A-D zoom into specific resonance structures to highlight the two-di m ensional lineshapes. 


frequencies, as well as by the detuning of the optical field. Experimentally this typically leads to 
relevant time scales on the order of hundreds of microsecond^. Moreover, the evolution can 
be slowed down artificially by adiabatically changing the ratio of trap frequencies. These time 
scales allow for the applications of pulses that, on the one hand, are short compared to the system 
dynamics, and, on the other hand, are long enough to spectrally resolve individual transitions. 

Whereas often the information contained in single quantum coherence signals is also available 
in technically less demanding Ramsey interferometric measurement sequenceP, this is not the 
case for the present simulation study. In a weak-field approximation, phase-cycling enables us to 
isolate those events where exactly N excitations are exchanged with the interaction pulses - of 
which no clear analogue exists in the context of a Ramsey experiment. 


B. Coherence decay of distinct quantum phases and evolution of the atomic 

insulator state 

We now present simulations of the three-pulse sequence ([^ with the level scheme in Fig. [T] for 
the ratio A/g = lOatA^ = 2. The initially prepared state |phSF) has strong overlap with the 
manifold ground state |(l/o )• Small corrections from the first single spin excitation submanifold 
become visible through oscillations during ti. After two interactions with the first pulse, the 
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strongest contribution to the generated state |atl) is given by the highest excited state |^^ 2 )- In 
the subsequent time interval t 2 , the decay of the coherence |atI)(phSF| can be observed in the 
two-dimensional spectrum. 

Fig. [^displays the two-dimensional spectrum = 0,^2, ^ 3 ) (where we omit the site in¬ 

dex j in ([SpO]), since the fluorescence signal is independent of j for N = 2). Our simulations were 
carried out using the qut ip packageP^. The parameters used in our simulation are ^ MHz, 
/3 = 5 kHz, A = 50 kHz, g = 5 kHz and reflect those of a recently reported experimental. Xq 
obtain realistic results, we further include the effect of the dominant error source - collective 
dephasing of the qubits due to fluctuations in the magnetic field that determines the common elec¬ 
tronic resonance frequenc3lAJ 34|35 |^ described by the dephasing operator af^ with decay constant 
7 = 0.5 kHz. Another relevant noise source is caused by fluctuations of the laser intensity of the 
driving field^, which we do not consider here. In the following, we analyze the position and 
two-dimensional lineshapes of cross-peaks, which are direct traces of the character (localized or 
delocalized) of the involved states. 

We identify four strong resonances in Fig. which are labelled A,B,C,D. All peaks can be 
uniquely attributed to a specific excitation pathway by extracting the corresponding resonance 
frequencies along O 2 and ^ 3 . Resonances A and C at r 22 ~ 672 kHz originate from the coher¬ 
ence 1 ^ 2 ) (^0 I during the f 2 -evolution. Hence, it is those two peaks that trace the evolution of 
the “macroscopic” coherent superposition of the two distinct phases (|phSF) + |atl))/v^. The 
stronger resonance A at Os = 0 signals the evolution of the population |^^ 2 )(^' 2 | during The 
states and ) are populated during due to the weak residual coupling of the spin-phonon 
interaction. Two less pronounced resonances at O 3 ~ 350 kHz reveal the coherent evolution of 
superposition states that involve the coherences |^' 2 )(^^f^^|- 

Furthermore, the two-dimensional lineshapes of the resonances can be used as a tool to identify 
the nature of the underlying states. The ability to analyze the two-dimensional lineshape^ rep¬ 
resents another distinct advantage of multi-dimensional spectroscopy, and becomes particularly 
useful in the present context. The influence of the collective dephasing process increases quadrat- 
ically with the difference of occupation numbers that contribute to the superpositiorP^H^, i.e., we 
have 0(f)|atI)(phSF| ~ and therefore is most pronounced for the decay of electronic co¬ 

herences between the ground state manifold and the highest excited state. When both states contain 
exactly the same number of spin excitations, they evolve in a decoherence-free subspacd^. Here, 
however, the non-zero coupling between spins and phonons perturbs the symmetry of the subspace 
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and will always lead to finite linewidths. Hence, both resonances A and C are broadened strongly 
along Q 2 (where the discrepancy of spin excitations is two), and much less along ^3 (where the 
states differ by no more than a single spin excitation). As discussed before, the resonance A 
represents the population |^ 2 )(^' 2 | during i.e., it evolves in a decoherence-free subspace. In 
contrast, the resonance C is also broadened along O 3 , but due to the smaller difference in the spin 
population between and 1 ^ 2 ). the broadening is less pronounced compared to the strong 

broadening along ^ 22 - 

Resonances B and D at ^2 ~ 350 kHz pertain to processes in which the two interactions with 
the first pulse only create a single spin excitation - again, due to the residual spin-phonon coupling 
- and the system therefore evolves in a coherence between one of the two single-spin states 
and ) and the initially prepared ground state )• resonance B, the second pulse then 
creates two spin excitations on the bra side of the density matrix, such that the system evolves 
in the coherence during t^. Consequently, the difference of spin excitations remains 

unchanged during both propagation intervals, and the resonances are broadened equally along both 
0,2 and O 3 . Finally, the resonances in the close-up D represent a population during t^, such that 
they are only broadened along ^2 - just like resonance A. 


C. Time evolution of populations of the two quantum phases 

We finally analyze the signal t 2 = 0, Os), Eq. (|^, for the case of negative detuning 

A = —50 kHz. Here, the two-spin state |5'2) forms the ground state, and the zero-spin manifold 
represents the highest-excited states [see Fig.[^)]. The signal correlates the frequencies of 
the two free evolution periods when the system is fully described by populations. Consequently, 
the main resonance can be found at zero frequency, with a broadening solely due to the finite 
propagation time. The signal measures the deviations from the \/^/g\ —> 00 limit, as the four 
resonances that can be found at nonzero frequencies all stem from deviations of the ground state 
|^' 2 ) from |atl), and of from |phSF). The resonances at Oi ±350 kHz yield the dominant 
correction to |atl), and the resonance pairs at Os ±600 kHz those to |phSF) in the simulated 
realization of the polaritonic Hamiltonian. 
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FIG. 4. Absolute value of the two-dimensional spectrum = 0,fi3)|, see Eq. (|^, of the 

pulse sequence ([^ in a two-ion chain, with the Hamiltonian ([^ and parameters = 1 MHz, /3 = 5 kHz, 
A = -50 kHz, 5 = 5kHz. 

IV. SUMMARY & OUTLOOK 

We have proposed a nonlinear measurement scheme to probe the nonequilibrium dynamics 
of polaritonic quasiparticles in trapped ion chains. Suitably chosen laser pulses can drive the 
system between an atomic insulator and a superfluid phase, and create superposition states of the 
two phases. On top of the identification of the contributing eigenstate frequencies, the analysis 
of the two-dimensional lineshapes in the measurement setup allows for a clear identification and 
characterization of the states, since the dominant noise source is sensitive to the number of spin 
excitations present in the system. Hence, the broadening of the resonance along specific axes 
indicates the discrepancy of electronic excitations between the two coherently superposed states. 
The variation of the third time delay, which was kept fixed in our present study, would further 
allow to monitor transport and relaxation processes within the A^-excitation submanifold. 

The measurement scheme proposed in this work does not require individual laser access to 
the ions and is therefore easily scalable towards larger N. Local information on the excitations 
can still be achieved by spatially resolved collection of the fluorescence light, which is standard 
practice in ion trap experiments. The distinct couplings of the ions to different collective modes 
are then reflected in the locally collected measurement signals. 

The system we studied - a collection of two-level atoms that couple to harmonic motional 
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modes - could also be interesting for experimental studies of the role of the vibrational backbone 
for exciton transport, which is believed to represent an important mechanism for biomolecular 
transport processe^^^^. To this end, one would model the direct exciton couplings, which is 
lacking in the polariton model, by adding effective spin-spin interactionP^l. in the context of 
biomolecules, an electron-phonon coupling that preserves the number of electronic excitations, 
i.e., a polaronic coupling, would be more appropriatd^ than the polaritonic coupling considered 
in the present study. 
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Appendix A: Derivation of the interaction operator 


In this appendix, we derive Eq. a. Our derivation closely follows standard approaches in 
the trapped-ion literature, which can be found, for instance in refs.^®^. Its starting point is the 
Hamiltonian 


N 


N 


N 


Ho{t) = '^uJkalak + ^ tu {a\ai^a\ak\ -Fa; 


k=l 


Q 


k,i=i 

k<i 


k=i 


N 


_ ^-i{kLXk-(^Lt+(l))^ ^ 


(Al) 


k=i 


which describes the interaction of the ion chain with the external laser field of frequency ujl, 
wavevector kL = which points along the x-axis, and phase (p. We assume spatially homoge¬ 
neous illumination of all ions, whose optical resonance frequency is given by uj. Note that due to 
the large separation of the ions (see section [TT aI ), the spatial dependence of the laser field needs to 
be taken into account. We first diagonalize the vibrational Hamiltonian, 


N 


N 


N 


^ ^ LO]^Qj^CL]^ + ^ ^ t]^i CLCL ^ ^ 


(A2) 


k=i 


k,i=i 

k<i 


n=l 


and express the motional dependence as a function of the Lamb-Dicke parameters r]n = kL^l/2mh'n 
oJMI 


exp 


{ikiXk) = exp i'^rjniCknbn + 4n^n) > 


(A3) 
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where the spatial orientation of the laser is adjusted such that only the vibrational modes along one 
direction are addressecP^. 

The excitation frequency is set on resonance with the red sideband corresponding to the vibra¬ 
tional eigenmode I'm, i-e., lol = io — Vm- Next we employ a rotating wave approximation and 
change to the interaction picture, employing the transformation 


U\it) = exp I —i 




0 , 


(A4) 


resulting in 


Hiit) = 


k 


exp ^ Vnicknbne-^^-^ + j + H.c. (A5) 


We shall consider the Lamb-Dicke limiP^, iin\J{ipn + bhy) <C 1, which is reached when the am¬ 
plitude of the ion’s motion is much less than 1/A:. In this case we may approximate the exponential 
in Eq. ( |A5| ) by its Taylor series to linear order, 

exp H r]nicknbne~"'"* + J 1 + z Vn{cknbne~"'"^ + . (A6) 

\ n / n 

Retaining only non-oscillatory terms in Eq. (|A5|), we then arrive at 


H,(t) - ) bm + H.c. 


N 


(A7) 


,fc=i 


For m = 1 we finally obtain Eq. by identifying Ckicr'^ and Q = rjiQ. 


1. The case N = 2 

In the present of case of a two-ion chain, we have two delocalized eigenmodes, which are given 
by the operators 


^2 — + ^2) ) 


(A8) 


the center-of-mass (high-energy) mode, and 


bi = (ai - a2), 


(A9) 


the breathing (low-energy) mode. Hence, for = 2 the coefficients in Eq. (|^ are cn = —C 21 = 
l/\/2. 
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